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ABSTRACT 

Aims. We study the nature of Gomez's Hamburger (IRAS 18059-3211; GoHam), a nebula that has been proposed to be a post-AGB 
object. Such a classification is not confirmed; instead, we argue that it will be a key object in the study of disks rotating around young 
stars. 

Methods. We present high resolution SMA maps of CO 7=2-1 in Gomez's Hamburger. The data are analyzed by means of a code 
that simulates the emission of a nebula showing a variety of physical conditions and kinematics. 

Results. Our observations clearly show that the CO emitting gas in Gomez's Hamburger forms a spectacular disk in keplerian rotation. 
Model calculations undoubtly confirm this result. The central (mainly stellar) mass is found to be high, ~ 4 M© for a distance of 500 
pc. The mass and (relatively low) luminosity of the source are, independent of the assumed distance, very diff'erent from those possible 
in evolved stars. This object is probably transitional between the pre-MS and MS phases, still showing interstellar material around the 
central star or stellar system. 
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1. Introduction 

Gomez's Hamburger (IRAS 18059-3211; hereafter GoHam) is a 
southern nebula with a very spectacular appearance in the optical 
(see |Ruiz et al. 11987) and HST images in press release number 
2002/19): an obscuring lane separates two flat, bright regions, 
presumably illuminated by a central star that remains hidden by 
the equatorial disk. The nebula was originally identified on a 
plate by A. Gomez, in 1985. The central star has been classified 
as an AO-type star, from spectroscopic analysis of the scattered 
light ( [Ruiz et aL]|1987| ). Its spectral energy distribution (SEP ) 
shows two maxima, in the optical and FIR ( [Ruiz et aL]|1987 ). 
Some emission excess at 2-10 jim also appears, probably due 
to relatively hot dust. Similar features are often detected in pre- 
MS and post-AGB sources, they have been proposed to be due 
to reservoirs of material relatively close to the star, i.e. rotat- 
ing disks (e.g.jCh iang & Goldreich| 19971 [Wood et al.|2002[[viT| 
Winckel 2003 1 ). In our case, the relatively low intensity of the op- 
tical maximum is due to strong absorption by grains in the disk, 
seen almost edge-on, the FIR maxima being due to reprocessed 

light 

|Ruiz et aL] ( [1987| proposed that GoHam is a post-AGB neb- 
ula, mainly based on similarities with some of such objects, like 
the Red Rectangle and M 1-92, and the absence of clear associ- 
ation with interstellar clouds. We also note the absence of con- 
spicuous emission components in the Balmer lines, which are 
typical of protostars. 
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We present CO /=2-l high-resolution mapping of GoHam 
that undoubtly indicates that the nebula is essentially a disk in 
keplerian rotation. As we will see, the rotation velocity yields 
a relatively high mass for the central star(s). We will argue that 
the central stellar mass and luminosity are not compatible with 
those acceptable for post-AGB sources, but are very similar to 
those expected in pre-MS A-type stars. 



2, Observations 

The observations were obtained with the Submillimeter Array 
(SMA) on the night of 2006 June 6 (UT), with the array in 
the extended configuration. All eight antennas were operational. 
The J2000 coordinates used for IRAS 18059-3211 were R.A. = 
18:09:13.37, Dec = -32:10:49.5. The receivers were tuned to 
place the ^^CO /=2-l transition at 230.53797 GHz in the upper 
sideband, and the ^^CO 2-1 transition at 220.39868 GHz in the 
lower sideband. During the track the precipitable water vapor 
above the array ranged from 1.2 to 2.0 mm. The DSB receiver 
temperatures ranged from 100->140 K near source transit, to 
160->190 K at the lower elevation limit of 15 degrees. The cor- 
relator was configured to give a spectral resolution of 203.125 
kHz (~0.26 kms"^) per channel in the regions of the IF con- 
taining the CO lines, and a coarser 3.25 MHz (-4.23 kms"^) 
per channel in the remaining regions of the full 2 GHz IF band- 
width. Ganymede was used as a flux calibrator, 3C279 was the 
bandpass calibrator, and observations of the quasars 1911-201 
(2.3 Jy) and 1924-292 (4.4 Jy) were made every 20 minutes, to 
track changes in the instrument's complex gain. IRAS 18059- 
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3211 was observed from 08:45 to 15:00 UT, and 4.6 hours of 
on-source data was acquired. Projected baseline lengths ranged 
from 12 to 140 kA. 

The MIR package was used for calibration. Antenna-based 
bandpass spectra were derived from the crosscorrelation data, 
and applied to flatten the spectral response. Next, pseudo- 
continuum data were derived by vector averaging the spectral 
channels in the region of the bandpass free of line emission. 
This pseudo-continuum channel was used to derive the fluxes 
of the quasars. The interleaved quasar observations were then 
used to calibrate the IRAS 18059-3211 data. The data were then 
exported to the AIPS package for mapping. Continuum emis- 
sion was clearly detected in our source, and it was subtracted 
from the channels showing line emission before mapping (see 
Fig. 1). The measured continuum flux at 1.3 mm is ~ 0.3 Jy; its 
deconvolved image is about V wide and elongated in the direc- 
tion of the nebula's major axis, at a position angle of about -5 
degrees (measured from north to east). The J2000 coordinates 
of the continuum source center are R.A. = 18:09:13.42, Dec = 
-32:10:50.0. 

From checks on the positions of the calibrators, we estimate 
that the accuracy in the absolute coordinates is ~ 0!'3. 

We present here the ^^CO /=2-l maps. Figs. 1 and 2. The Icr 
noise in the channel maps is 0.1 Jybeam"^ The half-maximum 
size of the clean beam is 1.47x1.09 arcsec, the RA. of the ma- 
jor axis being -10 degrees. The corresponding conversion factor 
from flux units to Rayleigh- Jeans-equivalent brightness temper- 
ature is 14.3 K/(Jy/beam). In Fig. 3 we show the line profiles of 
^^CO and ^^CO /=2-l resulting from integration over the emit- 
ting region. The ^^CO /=2-l observations, as well as ^^CO 3-2 
and 6-5 data obtained in a recent observing run, are still under 
analysis. 

We also show in Fig. 1, first panel, and image obtained from 
the HST archive. This image was taken on April 12, 2006; HST 
project 10603, P.I.: D. Padgett. The NICMOS camera was used, 
with the Fl lOW filter (wide filter centered on 1.10 microns) and 
an exposure time of 768 seconds. The position and orientation 
of the HST image is compatible with our mm- wave maps within 
the astrometric uncertainties. 

3. CO emission model 

We have analyzed our data using a model that simulates the 
emission of the nebula for a variety of physical conditions and 
kinematics; it was soon found that Gomez's Hamburger is well 
modelled by a disk in keplerian rotation whose width increases 
with the distance to the center. The code and general features of 
the model are very similar to those described by [Bujarrabal et al. 
( |2Q03l[2005] ). 

In the line emission model, the populations of the CO rota- 
tional levels are given by a single excitation temperature, which 
is assumed to be equal to the kinetic temperature. As discussed in 
our previous works, these approximations are valid for the typi- 
cal densities actually found from our model fitting (> 10^ cm"^) 
and the high opacities expected for the ^^CO lines (from com- 
parison of the data presented here with our preliminary ^^CO 
maps, see Fig. 3). We assume a shape of the CO cloud and a 
spatial distribution of the velocity, temperature, density, and CO 
relative abundance. Given these parameters, our code calculates 
the brightness for a number of lines of sight solving the full ra- 
diative transfer equation. Opacity eff'ects and velocity shifts are 
accurately taken into account. Such a brightness distribution is 
convolved with the cleaned telescope beam, and images with the 
same units as the observed ones are produced. 



The continuum emission at A = 1.3 mm is treated follow- 
ing 'Bi^OTabaretaL (2005 ); note that its contribution is not sub- 
tracted in the observed velocity-position diagram. Fig. 2. We as- 
sume the detected mm-wave continuum to be due to optically 
thin dust emission, since the measured intensities are compati- 
ble with the extrapolation of the FIR dust emission ( [Ruiz et aL 
1987|. The free-free emission in young A-type stars (Skinner 
1993| ) is usually much weaker than the level detected in 



et al. 



our data; in the Red Rectangle, a post-AGB source that could 
have similar properties as GoHam, the free-free contribution to 
the continuum at this frequency is also much lower than that of 
dust emission (e.g. |Men'shchikov et aL]|2 002 ). In any case, the 
treatment of the continuum is the same if there is some signif- 
icant contribution from optically thin free-free emission. Note 
that the very probably low optical depth of the dust continuum 
at this frequency simplifies the treatment and avoids discussions 
on the dust characteristic temperature. The spatial distribution 
of the dust emission coefl&cient is described in the code by a 
three-dimension gaussian function around the nebula center; the 
corresponding parameters are estimated from normalization of 
the resulting total flux and extent with respect to the data (Sect. 
2). The efl'ects of the dust continuum emission are however very 
small and hardly visible in the line maps. 

The disk shape is defined by its inner and outer radii (7?in, 
Rout) and thickness (H) at these radii. 

The macroscopic velocity is simply a keplerian rotation, de- 
fined by the tangential velocity at a given distance; the data 
are not compatible with the presence of significant expansion 
or infall velocities. The local velocity dispersion is assumed to 
be composed of the thermal dispersion (given by the kinetic 
temperature distribution, see below) and turbulent movements. 
However, including turbulence was not found to be necessary in 
our case to fit the data, whose low velocity dispersion indicates 
turbulent velocities < 0.1 kms"^ 

The total density n and temperature T of the disk gas are 
assumed to vary with potential laws. Thus, in our model: n(r) 
= n{RoWolrr^ and T{r) = T{Ro){RolrT\ where n{Ro), T(Ro), 
an, and ar are free parameters. 

We assume a CO relative abundance X(CO) = 10""^, usually 
found in interstellar condensations, but somewhat smaller than 
the values typical of young PNe. Note that, due to the assump- 
tion of a thermalized level population, changes in the abundance 
and density are compensated for exactly in the model predictions 
when the product nX remains constant. 

To account for the increase in intensity observed to come 
from the edges of the disk (Fig. 1), we also included possi- 
ble variations of the density and temperature in them. Since the 
^^CO lines are optically thick and this limb brightening does not 
appear in our ^^CO data, such a phenomenon must be mostly 
due to a higher temperature and a lower density in the edges 
(which is in any case expected in this kind of disk). We have 
successfully fitted the data assuming that the outer 20% of the 
cone (in the axial direction) shows a temperature increase by a 
factor 2 and a decrease in density by a factor 10 (with respect to 
the general potential laws). We note that the brightness increase 
is less remarkable in the east part of the disk than towards the 
west due to significant self absorption. Some limb brightening is 
also found in the outer edge of the cone (in the disk plane), then 
implying variations of the same kind in these regions: our best 
fit includes an increase of the temperature, by 20%, in the outer 
very diff'use part of the disk (outer 25%), in which the density 
already shows a decrease by a factor of 10. We also assumed 
an increase of the temperature of 50% in the innermost regions, 
within 7?in and 27?in. Although we are not comparing here the 
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Fig. 1. Channel maps of the ^^CO J=2-l line from Gomez's Hamburger, continuum has been subtracted. First contour and contour 
step are 0.2 Jy/beam; negative values are represented by dashed contours. The LSR velocity in km s"^ for the center of each channel 
is indicated in the upper left corner. The J2000 coordinates for the reference position, the cross in the maps, are R.A. = 18:09:13.37, 
Dec = -32:10:49.5. The first panel shows an HST/NICMOS image of this source. In the second panel we represent the map of the 
velocity-averaged CO intensity (first contour and step: 0.05 Jy/beam). The last panel gives the continuum intensity distribution (first 
contour and step: 0.03 Jy/beam) and the clean beam half-intensity size (black ellipse). 



physical conditions we find from our simplified fitting with disk 
models (to be done in a forthcoming paper), we note that such 
an increase of the temperature in disk edges is expected from 
theoretical results (e.g. Dulle mond et al.|2()07] ). 

Note that the observed disk emission is not exactly symmet- 
rical in velocity and position with respect to the rotation axis 
(RA. ~ 85°; see Figs. 1, 2); our model assumes axial symmetry 
and only can reproduce an average of both parts. The emission 
from the outer, very diff'use edge of the disk (spatial off'sets ~ ± 
5'', shift relative to systemic velocity ~ ± 1.2 km s"\ see Fig. 2), 
in particular, presents a noticeable asymmetry between the north 
and south parts of the disk. The asymmetry found in our maps 
with respect to disk equator (P. A. 5°) is due to the slight in- 
clination of the nebula axis with respect to the plane of the sky 
and well reproduced by our model. 



The fitted values of the parameters describing the structure, 
kinematics and physical conditions in the model nebula are given 
in Table 1. The parameter values are scaled with the distance; 
we take Z) = 500 pc as standard value following our final con- 
clusions on the possible values of the stellar mass and total lumi- 
nosity (Sect. 5). The density ranges between 1.5 10^ cm"^ (close 
to the central star) and 1.5 10^ cm"^, in the outer regions (ten 
times smaller in the very diff'use outermost edge). The temper- 
ature ranges between 65 K (close to the star) and 20 K (in the 
coolest equatorial regions). 

We note that some of the fitted parameters are almost di- 
rectly given by the observational features, their values are there- 
fore quite accurate and not model dependent. The model fitting 
of course confirms the values that we directly measure from the 
maps. The velocity variation with the distance to the star is ob- 
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Table 1. Structure and physical conditions in the molecular disk in Gomez's Hamburger, derived from our model fitting of the ^^CO 
J=2-l data. Dependence on the assumed distance is given in the relevant cases. Other parameters of the modeling are also given. 



Parameter 


Law 


Values 


comments 


Outer radius 


/?out = 4.75 1016 (^) cm 




Inner radius 


/?in = 3.5 10^5 (^) cm 




Disk 
thickness 


linear 


H{R,,,)^23 1016 (^) cm 
H{Rtn) = 5A 10^5 (^) cm 




Tangential 
velocity 


Vtocl/V^ 
(keplerian) 


yt(/?out/2) = 1.6kms-i 
central mass: 4 (^^) M© 




Temperature 


T oc 1/r^^ 


r(i?out/2) = 21 K 
ar = 0.4 


plus increase in disk edges 
(Sect. 3) 


Gas density 


noc l/r^'^ 


n(R,J2) = 2.5 106(^)cm-3 

an = 1 


plus decrease in outer edges 
(Sect. 3) 



Other parameters 


Law Values 


comments 


Axis inclination from the plane of the sky 


5° 


from optical and CO data 


Axis inclination in the plane of the sky (PA) 


85° 


from optical and CO data 


Distance 


500 pc 


from possible values of the luminosity and mass (Sects. 4 and 5) 


CO relative abundance 


constant 10""^ 


from various CO data (Sect. 3) 


LSR systemic velocity 


2.7kms-i 


from CO data 



viously kepler-like and the measurement of the central mass is 
accurate, see Fig. 2 and further discussion in Sect. 4. The diam- 
eter of the disk is also given by the observations, about ~ 13'', 
equivalent to ~ 10^^ cm at Z) = 500 pc. The total width of the 
disk (at least its maximum value, at the end of the disk) can also 
be measured from our maps, we find a maximum width ~ 3'\ 
2.3 10^^ cm at 500 pc. The data also indicate that the disk width 
increases more or less linearly with the distance to the center, the 
result being a flaring or quasi-flaring disk. Although such a thick 
disk cannot be considered to be exactly 'keplerian', the variation 
of the tangential velocity does follow the law. 

Since the ^^CO /=2-l line is found to be optically thick, the 
dependence of its intensity on the density and CO abundance is 
small, and these parameters cannot be well determined from our 
data. The values used here have been chosen to be compatible 
with our preliminary ^^CO maps, though no accurate fitting of 



our observations of both lines is attempted in this paper. In par- 
ticular, we cannot give here an accurate measurement of the disk 
total mass. On the other hand, the fact that the line is optically 
thick allows us to determine the temperature that, regardless of 
selfabsorption eff'ects, must be similar to the measured bright- 
ness temperature. 

The inner radius of the CO-emitting disk is also not well 
determined from our data. We detect CO emission up to ± 3.5 
km s"^ from the systemic velocity. For the assumed velocity law, 
this means that we detect emission from regions farther than ~ 
0'.'5. This value is similar to the beam radius, as expected, be- 
cause much smaller regions are strongly diluted within our spa- 
tial resolution and more difficult to detect. On the other hand, we 
do not see a clear central minimum in our maps, confirming that 
the inner disk radius cannot be larger than the beam. Regions 
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closer to the center than ~ (7/5 may also be molecule rich, but 
they are not well probed by our data. 

Our synthetic maps essentially depend on the nebula struc- 
ture, velocity field and temperature distribution. In this work, we 
focus on the determination of the rotation velocity and the cen- 
tral (mainly stellar) mass, which involves in fact only one free 
parameter, Vt at a given point. The distance is widely discussed 
in Sects. 4 and 5. Accordingly, the model remains very simple. 
The determination of the physical conditions (density, tempera- 
ture, CO abundance) will be improved in a forthcoming paper, 
from fitting of several rotational transitions of ^^CO and ^^CO 
and a detailed discussion on the physical meaning of the results 
and their compatibility with theoretical considerations. 

4. The mass and luminosity of Gomez's Hamburger 

Our observations clearly show that the CO emitting region in 
Gomez's Hamburger is a disk, strikingly coincident with the 
equatorial dark lane and bright rims seen in the optical. The CO- 
rich disk is even slightly larger than the optical image and very 
probably represents the total nebular component of this source 
(Fig. 1). The position- velocity diagram along the equatorial di- 
rection clearly shows the butterfly-like shape characteristic of 
keplerian rotation (Fig. 2). Our model of emission from a ro- 
tating disk reproduces very closely the CO line observations of 
GoHam, see Fig. 2, confirming that the Gomez's Hamburger 
nebula is a disk in keplerian rotation. 

The optical image and the double-maximum SED of GoHam 
(Sect. l,|Ruiz et al.|1987)) strongly support the identification of a 
disk seen practically edge-on. The energy radiated by our source 
shows a double maximum in the optical and FIR, with also some 
excess at about 2-10 yum, which is characteristic of a disk rotat- 
ing around a relatively hot star. Both young and post-AGB stars 
surrounded by disks of dusty material show SEDs of this kind 
(e.g. [van Wincke l 2003 ; Wood et al.||2002| ), in which the short- 
A maximum comes from the stellar emission, absorbed and/or 
scattered by the grains, while the FIR maximum is due to the 
emission of the disk dust, heated by the stellar radiation. 

Our model fits the rotating disk structure for a central mass 
Ms(Mq) ~ 8 Mq, for a distance D = I kpc; Ms varies propor- 
tionally to the assumed distance. It is important to note that this 
result is very model-independent. Since the disk is seen almost 
edge-on, corrections to the tangential velocity and distance to 
the center due to the geometry are negligible. We can directly 
see from our maps and velocity-position diagram that gas at a 
distance of 6 arcseconds (equal to 9 10^^ cm if Z) = 1 kpc) ro- 
tates at 1.1 kms"^ Errors to these figures, directly taken from 
the maps, are smaller than ~ 10%, including the correction due 
to projection eff'ects (the correction to the observed velocity for 
an inclination of 10°, see Sect. 3, is of less than 2%). These val- 
ues of the velocity and distance also yield, assuming keplerian 
rotation, a central stellar mass M^(Mq) = SMq x Z)(kpc). The 
choice of other velocity/position pairs yield the same mass, be- 
cause, as previously mentioned, the velocity field is very exactly 
keplerian. Such a mass is very high for a post-AGB star, since 
the central stars of planetary nebulae are expected to have masses 
smaller than 1 M©, see e.g. Pottasch (1984 ), Blocker (1995 ). 

On the other hand, the total luminosity of GoHam, inte- 
grating all wavelengths, is moderate, ~ 35 Lq for Z) = 1 kpc, 
since the apparent bolometric luminosity is 1 . 1 10~^ erg cm"^ s" 



(from careful study of the SED by |Ruiz et al.|1987| ). [RuizetaLl 
(1987]) compared the (apparent) SED of GoHam with those of 
two well known post-AGB objects, the Red Rectangle and M 1- 
92 (see their Fig. 4). GoHam is more than ten times weaker at 



all wavelengths than M 1-92 and about 100 times weaker than 
the Red Rectangle. M 1-92 has an absolute luminosity of 10^ 
Lq and a distance of about 2.5 kpc; the Red Rectangle has a 
luminosity -3.5 10^ Lq, assuming a distance of 710 pc (e.g. 
[Bujarrabal et al.|2001t[Men'shchikov et al.|2002) . Both are typi- 
cal post-AGB stars, and other objects in this evolutionary phase 
show very similar absolute luminosities. We also note that the 
mass of the binary system in the center of the Red Rectangle is 
well known, 1.5 Mq, measured from the (almost) keplerian disk 
around it. 

Accordingly, for GoHam to have a post-AGB luminosity, at 
least, ~ 3 10^ L0, it should be placed at a surprising distance D 
-10 kpc. This would imply an exceptional central mass Ms ~ 
80 Mq, as well as a distance to the galactic plane ~ 1 kpc, par- 
ticularly unexpected for such a massive object. Even assuming 
that there is a binary in the center, we should need something 
like two 40Mo stars. 

If we assume that the stellar component is a binary and most 
of the mass is in the less evolved star, it should have a mass of 
almost 80 Mq. Then, the evolved star should have had an initial 
mass still larger (since it evolved faster), leading to still more 
surprising conclusions. 

These conditions can be somewhat relaxed if we take into 
account that the stellar sho rt-/l emiss ion is absorbed by dust be- 
fore it leaves the nebula. R uiz et al. (^1987 ) suggest a reddening 
E(B - V) ~ 0.56; we must take into account that the interstel- 
lar e xtinction measured in the direction of our source is ~ 0.5 
mag ( Dobashi et al.|2005 ). In any case we note that a significant 
absorption in the optical in the direction in the equatorial plane 
is clearly suggested by the nebula image. This eff'ect would im- 
ply a strongly anisotropic emission in the optical and UV, being 
minimum in the direction of the equator (i.e. toward us). We can 
easily check that the correction due to the above extinction value 
is of the order of a factor ten (at short A). Larger correction val- 
ues cannot be obtained for reasonable geometrical distributions 
of the circumstellar dust. Taking into account that the total emis- 
sion from the stellar component of the actually observed SED 
is about 10 times smaller than the dust emission component, we 
conclude that the correction to the total luminosity due to extinc- 
tion would increase it by a about factor two, up to 70 Lq at Z) = 1 
kpc. The 'post-AGB distance' estimated for this value of the lu- 
minosity would be ~ 7 kpc, and the central mass should be close 
to 60 Mq, still far too large. The correction of the extinction 
eff'ects is uncertain, but in fact can hardly be much larger than 
this value, since the fraction of the total stellar light absorbed by 
the dust lane and reemitted in the FIR cannot be negligible. For 
the disk aperture deduced from our model fitting, about 1/3 of 
the total stellar light is intercepted by the dust lane. Therefore, 
the total luminosity would be at most three times larger than the 
FIR luminosity (about 100 Lq, at 1 kpc), corresponding to the 
extreme case of negligible contribution of direct and scattered 
stellar light to the detected fluxes. Even in this case, the efl'ects 
on the mass are very moderate, and a too large central stellar 
mass of ~ 45 Mq would be required. A similar correction is 
found in the models developed by |Men'shchikov et"aL ( |2002 ) for 
the Red Rectangle, a post-AGB star showing a massive rotating 
disk that is seen almost edge-on. We note that in other post-AGB 
stars showing rotating disks, the extinction due to dust is always 
moderate and the eff'ects of the viewing angle on the luminos 
ity determinati on are expected to be almost negligible (e.g. 
|Winckel|2003l ). 

Even if we are now checking if GoHam can be a post-AGB 
object, we will also consider the efl'ects of the viewing angle 
on the luminosity determination predicted by models of disks 



van 
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around pre-main- sequence stars. Moderate effects not exceeding 
changes by a factor two in the total luminosity are found by fit- 
ting of the SEP using t he models for youn g stellar objects by 
Robitaille et al.| ( |2QQ7| ). [Wood et al.|p502t published a grid of 
calculations for different disk models. In extreme cases, these 
calculations yield a correction to the total luminosity of about a 
factor 10. Even adopting such a large correction, we would de- 
rive a too high value of the central stellar mass, ~ 25 Mq. We re- 
call that these models are not well adapted to our object, mainly 
if we assume by the moment that it is in the post-AGB phase. 
Finally we note that the reasons why even very large corrections 
to the luminosity estimate cannot yield mass values comparable 
to that of an evolved star are: a) the luminosity depends on Z)^, 
while the mass derived from the velocity pattern varies like Z), so 
large variations of the luminosity yield much smaller variations 
of the mass, and b) the mass of post-AGB stars are expected to 
be significantly smaller than 1 Mq, typically ~ 0.5 Mq, really 
very far from the values we derive. 

Another option is to assume that the central mass is not 
stellar, but corresponds to a very massive and compact nebular 
component. Young planetary nebulae are much less massive, at 
most ~ 1 Mq, and not so compact (see Bujarrabal et al. 200 1|. 
There are two yellow hypergiants, IRC+10 420 and AFGL 2343, 
surrounded by nebulae with about 5 Mq ( [Castro-Carrizo et aL 



[2007 ), but these nebulae are spherical and very extended, both 
objects are placed at less than 500 pc from the galactic plane, 
and both have luminosities ~ 5 10^ Lq. 

We note, moreover, that at a distance of ~ 10 kpc, the rotat- 
ing disk would be a monster, with a total radius of ~ 10^^ cm 
and a total nebular mass larger than 10 Mq, nothing comparable 
to known disks rotating around post -AGB stars (see ,BujarrabaI[ 
let al.|2005| [20071 Ivan Winckel|2003D . 



5. Gomez's Hamburger: A rotating nebula around a 
young pre-MS star 

We have seen that it seems impossible to reconcile our knowl- 
edge on the mass and luminosity of Gomez's Hamburger 
(GoHam) with the assumption that it is a post-AGB object. The 
discrepancy is so big that even the highest corrections to the to- 
tal luminosity and distance are far from yielding an acceptable 
result. 

Young stars may also be surrounded by CO-rich rotating 
disks (e.g. putrey et al.|2007| ). GoHam presents properties that, 
in fact, seem closer to those of rotating disks around pre-MS 
A-type stars than to those of post-AGB nebulae. The CO line 
profile (Fig. 3) is relatively symmetric and shows two peaks sep- 
arated by about 2 k ms~^ This structure is very sim ilar to that 
found in AeBe stars ( jPent et al.|2005t|Raman et al.|2 006, etc), 
particularly in stars isolated from interstellar clouds. However, 
post-AGB nebulae sh ow in general very wid e profiles with high- 
velocity wings (e.g. Bujarrab al et al.|[200T] ), due to outflowing 
gas. Even in the Red Rectangle, the only post-AGB object in 
which a keplerian disk has been well identified, and in 89 Her, 
the best post-AGB candidate to harbor such a disk after the Red 
Rectangle, the CO profiles are different, showing a narrow cen- 
tral peak plus wide wings. In these cases, the wings are due to 
expanding outer gas in the nebulae. We recall that no trace of 
expanding gas has been found in GoHam. 

We also note that the beautiful optical image of 
Gomez's Hamburger (from which its name comes) shows a flat, 
three-component structure re markably similar to optical images 
of disks around young stars ( [Grosso et al. 2003 [ Perrin et al. 
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Fig. 2. Velocity-position diagram of the ^^CO /=2-l line emis- 
sion from Gomez's Hamburger (upper panel) along a direction 
defined by RA.: -5°. Lower panel: model fitting of the observed 
velocity-position diagram. In both cases, the same contours as 
in Fig. 1 are used, but continuum is not subtracted in these dia- 
grams. 



2006; "^ ^tson et al.|2007 etc). The optical image of our source 
is however very different from the extended, prolate images 
of post-AGB nebu la, including that of the Red Rectangle, e.g. 
ICohen et al.| ( |2004l ). 

If we assume that the star is a young pre-MS A-type star, 
some conflict still persists between the expected values of the 
stellar mass and luminosity and those derived from observations, 
but it is certainly much less severe. For a distance of 1 kpc, we 
have derived a central stellar mass ~ 8 Mq and a total luminosity 
of about 35 L©, or ~ 70 L© for a strong extinction in the line 
of sight. If we assume that there is a binary stellar system in 
the center, a pre-MS primary with about 4 Mq and ~ 10000 
K of surface temperature (value derived from spectroscopic and 
continuum data by nRuiz et al.| 19871 ) should have a luminosity of 
about 140 L©, from comparison with our knowledge on similar 
young objects (see compilation of data on other pre-MS sources 



van 



and evolutionary tr ack calculations in, e.g., |Pietu et al.|2003] 
|de n Ancker et a l.|1997 ). 

Measured and expected values of the mass and luminosity 
are closer if we assume a shorter distance. For D = 500 pc, a 
mass of the primary of about 2 M© is derived from our CO data 
(assuming a binary system). For such a distance, the total lumi- 
nosity of GoHam derived from its SED is ~ 18 L©. On the other 
hand, an A-type pre- MS star with 2 M© is expected to show 
luminosities ~ 15 L© ( Pietu et aLl|2003[ van den Ancker et al. 
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Fig. 3. ^^CO and ^^CO J=2-l line profiles, obtained from inte- 
grating in our maps over the emitting region. 
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1997| ), compatible with the measured value. We note that, even 
in the pre-MS scenario, it is in principle necessary to assume 
that there is a binary system in the center of GoHam to be fully 
in agreement with theoretical evolutionary tracks. Otherwise, we 
should have to assume a shorter distance leading to luminosities 
under 10 Lq, and young stars with such a low luminosity are not 
expected to reach a surface temperature of 10000 K. 

The stellar mass requirement is somewhat relaxed if we 
take into account that the central nebular mass is not negligi- 
ble. From our continuum measurement, which is probably due 
to dust emission (Sects. 2, 3), we can infer the total mass for the 
innermost region (within ~ T')- From the simple formulation 
described by Natta et al. (2000, see details on the assumed dust 

properties in that paper), we derive a mass of about 0.9 x (^^) 
Mq. This may particularly aff'ect the mass of the companion, 
and our stellar system would be formed by a AO-type star with 
~ 2 Mq plus a cooler companion (probably a T Tauri star) with 
about 1 Mq, always for a distance of 500 pc. We note that in 
the case of extreme absorption in the line of sight (expected for 
some heavy disks around pre-MS stars. Sect. 4), which yields a 
total luminosity of about 90 Lq at 500 pc, only one central star 
with about 3 Mq could explain both the total luminosity and the 
rotation curve of GoHam. The location of GoHam in the H-R 
diagram, assuming that it is a young star and taking into account 
these uncertainties in the luminosity, is shown in Fig. 4; we also 
represent well studied AeBe stars and evolutionary tracks from 
|van den Ancker et al.1(|1998[ ). 

As noted by I Ruiz et al.| ( 1987 ), Gomez's Hamburger is not 
associated with a conspicuous interstellar cloud. But it is within 
a region of the sky with a noticeable extinction of about -0.5 
mag ( [Dobashi et al.||20()5j ). We also note that many young A- 
type stars are not clearly associated with clouds (e.g. Dent et al. 
2005| ), particularly in the case of low-luminosity stars, for which 
the pre-MS evolution is thought to be relatively slower. The low 
luminosity, the absence of conspicuous emission in Balmer lines, 
and the loose relation with interstellar clouds strongly suggest 
that our source is a relatively evolved star, perhaps close to en- 
tering the main sequence. 

We do not think that the central star of GoHam is a young 
main sequence star, similar to the well known debris-disk ob- 
jects. The absence of emission in the Balmer lines would support 



Fig. 4. Location of Gomez's Hamburger in the H-R diagram, 
thick vertical line, taking into account the uncertainty in the lu- 
minosity due to the viewing angle eff'ects. We assume a distance 
of 500 pc, see Sect. 5; the stellar temperature, ~ 10000 K, is rel- 
atively well known (see Ruiz et al. 1987 and Sect. 5). GoHam is 
compared with well studied AeBe stars and evolutionary tracks, 
from [van den Ancker et aL] ( |1998| ) . 



this possibility (although there is a population of A-type stars as- 
sumed t o be very young that do not show such emission features, 
see e.g. |The et al.|1994[|Manoj et al.| 20°06). But the disks around 
young MS stars are thought to disappear relatively quickly (i.e. 
in ~ 10^ yr) and the debris disks are not massive, yielding a mod- 
erate FIR excess and being very rarely detected in CO emission 
(e.g. Thi et a l. 2001 ; De nt et al.|2005] ). 

In summary, we conclude that GoHam is very likely a young 
star at a distance of about 500 pc, still showing signs of inter- 
stellar material around it. GoHam is probably a transitional ob- 
ject between the pre-MS and MS phases, similar to stars like 
HD 141569, still surrounded by relatively massive circumstellar 
disks (e.g. |Merm et al.|20"Q4{ Augereau et al. 1999). In any case, 
the gas disk orbiting the central star(s) of GoHam is exceptional: 
it is very massive and extended, and its study is favored by its 
observing angle (the disk is seen almost edge-on) and the lack 
of contamination due to surrounding interstellar clouds. We are 
accordingly convinced that GoHam will represent a milestone in 
the study of disks rotating around young stars. 
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